Spin Fluctuation Induced Superconductivity Controlled by Orbital Fluctuation 
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A microscopic Hamiltonian reflecting the correct symmetry of /-orbitals is proposed to discuss su- 
perconductivity in heavy fermion systems. In the orbitally degenerate region in which not only spin 
fluctuations but also orbital fluctuations develop considerably, cancellation between spin and orbital 
fluctuations destabilizes d I .2_ !/ 2-wave superconductivity. Entering the non-degenerate region by in- 
creasing the crystalline electric field, d x 2_ y 2-wscve superconductivity mediated by antiferromagnetic 
spin fluctuations emerges out of the suppression of orbital fluctuations. We argue that the present 
scenario can be applied to recently discovered superconductors CeTIns (T=Ir, Rh, and Co). 



Unconventional superconductivity has been one of cen- 
tral issues in the research field of strongly correlated elec- 
tron systems. Especially since the discovery of high tem- 
perature superconductivity in cuprates, much effort has 
been focussed on elucidating the mechanism of uncon- 
ventional superconductivity, clarifying that a crucial role 
is played by antiferromagnetic (AF) "spin fluctuations" 
(SF) Q. The importance of AFSF is widely recognized, 
for instance, in 4/- and 5/-electron superconducting ma- 
terials such as CeCu2Si2 and UPd2Al3 Q as well as 
in organic superconductors such as k- (BEDT-TTF) g. 
Thus, it is widely believed that a broad class of uncon- 
ventional superconductors originates from SF. 

In d- and /-electron systems, however, the potential 
importance of orbital degrees of freedom has recently 
been discussed intensively. In fact, "orbital ordering" is 
found to be a key issue for understanding microscopic as- 
pects of the charge-ordered phase in colossal magnetore- 
sistive manganites || . This orbital ordering is primarily 
relevant to the insulating phase, while "orbital fluctua- 
tions" (OF) should be significant in the metallic phase. 
Recently, the effects of OF have attracted attention, since 
it is hoped to provide a new scenario for superconductiv- 
ity H. Especially from a conceptual viewpoint, it is im- 
portant to clarify how superconductivity emerges when 
both SF and OF play active roles. 

As a typical material for investigating superconductiv- 
ity in a system with both SF and OF, let us introduce 
the recently discovered heavy-fermion superconductors 
CeTIn 5 (T=Ir, Rh, and Co) § with the HoCoGa 5 -type 
tetragonal crystal structure. Due to this structure and 
to strong correlation effects, the similarity with cuprates 
has been emphasized. In particular, AFSF also plays an 
essential role in the Ce-115 system, since it exhibits quasi 
two-dimensional Fermi surfaces || and the AF phase ex- 
ists next to the superconducting state in the phase dia- 
gram of CeRhi_ x Ir x In 5 [|J. In fact, a line node in the gap 
function has been observed in CeTIns by various exper- 
imental techniques [fl0|| . It may be true that supercon- 
ductivity itself is due to AFSF, but an important role for 
the OF has been overlooked in spite of the fact that both 
SF and OF are originally included in the ground-state 



multiplet of Ce 3+ ion. In actual materials, superconduc- 
tivity occurs in a situation where OF is suppressed, since 
orbital degeneracy is lifted by the effect of the crystalline 
electric field (CEF). Thus, we envision a scenario where 
OF controls the stability of the superconductivity even 
though it originates from AFSF. 

In this Letter, we investigate superconductivity based 
on the orbitally degenerate Hubbard model constructed 
by the tight-binding method pd| ] . Solving the gap equa- 
tion with the pairing interaction evaluated using the ran- 
dom phase approximation (RPA), we obtain several su- 
perconducting phases around the spin and orbital or- 
dered phases. In the orbitally degenerate region where 
both SF and OF are developed, it is found that singlet 
superconductivity is suppressed due to the competition 
between them, while triplet superconductivity is favored 
since they are cooperative in this case. When a level split- 
ting is included to lift the orbital degeneracy, dj^^-wave 
superconductivity due to AFSF is stabilized in the vicin- 
ity of the AF phase. Thus, we claim that AFSF-induced 
supercondutivity in Ce-115 systems is substantiated in 
consequence of the suppression of OF. 

In order to construct the microscopic model for /- 
electron systems, let us start our discussion from a lo- 
cal basis of the Ce 3+ ion. Among the 14-fold degener- 
ate 4/-electron states, due to the effect of strong spin- 
orbit coupling, only the j=5/2 sextuplet effectively con- 
tributes to the low-energy excitations (j is total angular 
momentum). This sextuplet is further split into a 
doublet and a Tg quadruplet due to the effect o f cub ic 
CEF, where the eigen states are given by |T7±)= ^/1/6|± 

5/2)-V576| T 3/2), \r$£)= ±S/2)+jlJ&\ T 3/2), 

(2) 

and |rg_^)=| ± 1/2). Here + and — in the subscripts 
denote "pseudo-spin" up and down, respectively, within 
each Kramers doublet. 

Now we discuss the relative positions of the energy lev- 
els of Tt and rg T ' while taking account of certain features 
of CeTIns. Since CeTIns has a tetragonal crystal struc- 
ture and a quasi two-dimensional Fermi surface || , it is 
natural to consider a two-dimensional square lattice com- 
posed of Ce ions. Due to the effect of anions surrounding 
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the Ce ion, it is deduced that the energy level of T7 be- 
comes higher than those of Tg's. Thus, in the following, 
the T7 orbital is neglected for simplicity. Note also that 
we need to include the effect of the tetragonal CEF, which 
lifts the degeneracy of the T^ T) . Although a mixing be- 
tween Tf and Tg 1 ^ generally occurs under a tetragonal 
CEF, such a mixing is expected to be small since the 
orbital has higher energy. Thus, in this situation, 
the effect of a tetragonal CEF can be included in terms 
of a level splitting e between the and orbitals. 
In order to explain the magnetic anisotropy observed in 
experiments, must be lower than T g 2 \ i.e., £>0 in 
CeTIns. The magnitude of e for each CeTIns compound 
will be discussed later. It should be noted that the above 
level scheme is consistent with the following two facts: 

(i) Experimental results for CeTIns exhibit a larger uni- 
form susceptibility for magnetic field perpendicular to the 
Celn 3 plane than for the parallel case j7|Jl^] . This signif- 
icant anisotropy is well explained under the assumption 
that Ty is not the lowest-energy state and e is positive. 

(ii) The band-structure calculation results suggest that 
the almost flat band corresponding to T7 appears above 
the Fermi level Jll[ . 

In order to include the itinerant features of the 4/- 
electrons, a simple way is to include nearest-neighbor 
hoppin g for the /-electrons by the tight-binding method 
p],^3[ . Although the hybridization with the In 5p elec- 
tronic states may be important, here such an effect is 
considered as renormalization of the effective hopping 
amplitude of /-quasiparticles. Further, by adding the 
on-site Coulomb interaction terms among the /-electons, 
the Hamiltonian becomes 

iarr'a i 

+ U^n iT in iTl + U'^n ila n i2 a', (1) 

ir icrcr' 

where /i r(T is the annihilation operator for an /-electron 
with pseudo-spin a in the orbital Ig at site i, a 
is the vector connecting nearest-neighbor sites, and 
ni T a=fi Tcr fiTa- The first term represents the nearest- 
neighbor hopping of /-electrons with the amplitude t* T , 

between Tg T ^ and T 8 T ' along the a-direction, given by 
til=— -s/3t*2= —VSt* ! =3*22=1 for a=x and t\ 1 =^/3t\ 2 = 
v3*2i = ^22 = l f° r a= y> respectively, in energy units 
where 4^=1. Note the positive sign of the first term in 
H, since the M-point, not the T-point, is at the bottom of 
the bands forming the Fermi surfaces (ll]] . Note also that 
the present t$ T , is just the same as that of the e g electrons 
H , but this point will be discussed elsewhere. The second 
term denotes the tetragonal CEF. In the third and fourth 
terms, U and U' are the intra- and inter-orbital Coulomb 
interations, respectively. In reality, U—U', since they 
originate from the same Coulomb interactions among /- 
orbitals in the j=5/2 multiplet, but in this paper, we also 
treat the case where Uj^U' in order to analyze the roles 



of SF and OF. Since we consider quarter-filling (one /- 
elecron per site), the model in the limit of e=oo reduce 
to the half-filled, single-orbital Hubbard model. 

Now, in order to investigate superconductivity around 
the spin and/or orbital ordered phases, we calculate spin 
and orbital susceptibilities, X s (q) and x°(q), respectively. 
Within the RPA, these are given in a matrix form as 

X B (q) = [i-C> s x(q)]- 1 x(q) ) (2) 
X°(q) = [i + ?7 x(q)]- 1 x(q), (3) 

where labels of row and column in the matrix appear in 
the order 11, 22, 12, and 21, these being pairs of orbital 
indices 1 and 2. Note that 1 is the 4x4 unit matrix. 
l> is given by Uf hn =UI 2t22 =U , ^12,12=^21,21=^'. oth- 
erwise zero, while U° is expressed as n=U 22 22 =U, 
U? h22 =U^ n =2U', U? 2il2 =U$ li21 =-U', otherwise zero. 
The matrix elements of x(q) are defined by X^api^) — 
-TEk, n Ga/l(k+q, iw n )G^ (k, iu> n ), where T is temper- 
ature, G^v (k, iw n ) is the non-interacting Green's func- 
tion for /-electrons with momentum k propagating be- 
tween fi- and (/-orbitals, and o;„=7rT(2n+l) with in- 
teger n. The instabilities for the spin- and orbital- 
ordered phases are determined by the conditions of 
det[l - C/ s x(q)]=0 and dct[l + C/°x(q)]=0, respectively. 

By using x s (q) an d X°(q), the superconducting gap 
equation is given as 

A«(k)=£v%-k')0(k')A«(k'), (4) 

k' 

where A«(k) =[Af x (k), A« 22 (k), A? 2 (k), Af x (k)]* is the 
gap function in the vector representation for a singlet 
(£ = S) or triplet (£ = T) pairing state, and the matrix 
elements of the singlet- and triplet-pairing potentials, re- 
spectively, are given by 

V^M = [(-3/2)^ s (q) + (l/2)VF°(q) + U s ] a ^, (5) 
V^M) = [(V2)^ s (q) + (l/2)^°(q) - U%^ p . (6) 

Here, the spin and orbital susceptibilities are included as 
W s {(i)=U s + U s x s (<l)U s and VF°(q)=-[> + J7°x°(q)t> . 
The elements of the pair correlation function 0(k) are 

given by </> Q/3iA1 „(k)=T J2 n Ga/i(k, -iu> n )G^ (-k, iw n ). 
The superconducting transition is obtained for each re- 
spective irreducible representation by solving Eq. (Q), 
where the maximum eigenvalue becomes unity. 

Here we mention the essential symmetries of the 
Cooper pairs in multi-orbital systems. The pairing states 
are classified into four types owing to spin SU(2) sym- 
metry and space inversion symmetry; (1) spin-singlet 
and orbital-symmetric with even-parity, (2) spin-triplet 
and orbital-symmetric with odd-parity, (3) spin-singlet 
and orbital-antisymmetric with odd-parity, and (4) spin- 
triplet and orbital-antisymmetric with even-parity. Note 
that SU (2) symmetry does not exist in orbital space, 
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since the lattice and the local wavefunctions rotate si- 
multaneously. For (1) and (2), /-electrons in the same 
band form the Cooper pair, while for (3) and (4), such 
a pair is formed only between different bands. Except 
for the special case in which two Fermi surfaces connect 
with each other, orbital-antisymmetric pairing is unsta- 
ble due to depairing effects destroy inter-band pairs. In 
fact, even after careful calculations, we do not find any 
region for (3) and (4) in the parameter space considered 
here. Thus, in the following, we discuss only the orbital- 
symmetric pair. 
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FIG. 1. (a) Spin and (b) orbital susceptibilities in q space 
for e=0 and U=U'=Q.9U m . (c) and (d) are for e—4. Insets 
indicate the Fermi-surface lines. 

Let us discuss first the effect of e on the correlation 
between SF and OF. Since we are considering a realistic 
situation corresponding to CeTIns, we restrict ourselves 
to the case of U=U' for a while. In Fig. 1, the principal 
components of x 8 (q) and x°(q) are shown in q space for 
U=U'=0.9U m , where U m denotes spin instability deter- 
mined from det[l — U s x(q)]—0. In the orbital-degenerate 
case with £=0, as shown in Figs. 1(a) and (b), the over- 
all magnitude of OF is comparable with that of SF. As 
for the q dependence, common structures are found in 
X s (q) and x°( c l)- Namely, there are two peaks around 
(7r/2,7r/2) and (7r/2,0) owing to the nesting properties 
of the Fermi surface (see the inset). Thus, the orbital- 
degenerate region is characterzied by competition be- 
tween SF and OF. With increasing e, the Fermi surface 
approaches a shape having the nesting vector (7r, 7t) (see 
the inset). For e=4, as shown in Fig. 1(c) and (d), OF is 
almost completely suppressed, and SF around (it, it) be- 
come dominant. Note here that the increase of e makes 
the lower energy state favorable and suppresses excita- 
tions to the upper energy state, indicating the suppres- 
sion of OF. Thus, even at this stage, it is understood that 
the suppression of OF leads to the development of SF for 
the paramagnetic system which is wavering between spin- 



and orbital-ordering. 

Next, we consider how a superconducting phase 
emerges when e is increased. In Fig. 2(a), the maximum 
eigenvalues for several irreducible representations are de- 
picted as a function of e for U=U'=2.5. The calculations 
are carried out for a fixed T=0.02, and the first Brillouin 
zone is divided into 128 x 128 meshes. In the orbital- 
degenerate region, several eigenvalues are very close to 
each other owing to multi-peak structures in x s (q) an d 
X°(q). With increasing e, as is easily understood from 
the growth of AFSF mentioned above, the eigenvalue for 
Bi g symmetry becomes dominant and finally at e~3, the 
Big superconducting phase is stabilized. On further in- 
creasing e, the AF instability eventually occurs, since the 
system asymptotically approaches the half-filled single- 
orbital Hubbard model. We emphasize that the increase 
of e brings transitions successively in the order of PM, 
Bi g superconducting, and AF phases. 
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FIG. 2. (a) Maximum eigenvalue vs. e for the respective 
irreducible representation at U=U'=2.5. (b) Schematic plot 
c/a vs. e to illustrate a comparison between our theory and 
actual CeTIns compounds. 

Based on the above calculated results, let us try 
to explain the differences among three Ce-115 com- 
pounds, CcRhlns (Neel temperature Tn=3.8K), Cclrlns 
(superconducting transition temperature T C =0.4K), and 
CeCoIn 5 (T C =2.3K) @. One property which distin- 
guishes these compounds is two-dimensionality, as ex- 
pressed by the ratio c/a between lattice constants. It 
is naturally expected that two-dimensionality becomes 
stronger in the order of Celrlxis , CeRhlns , and CeCoIn5 
JtJ. Since the magnetic state should be stabilized with 
increasing three-dimensionality, Celrlns should be most 
favorable for the occurence of antiferromagnetism among 
three compounds, but this is obviously inconsistent with 
experimental results. This inconsistency is resolved by 
introducing another important ingredient e. One can 
show that the increase of magnetic anisotropy just corre- 
sponds to the increase of e. Analysis of experimental re- 
sults for the anisotropy of magnetic susceptibilities [0JlJ] 
leads to the conclusion that e becomes larger in the order 
of CeCoIns, Celrlns, and CeRhl^. Taking into account 
the effect of e in addition to dimensionality, we arrive 
at the picture schematically shown in Fig. 2(b). Due 
to the enhancement of AFSF induced by increasing e, 
as shown in Fig. 1, it may be understood that CeRhIii5 
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rather than Celrlri5 is more favorable to antiferromag- 
netism. Namely, it is considered that CeRhlns with the 
largest e is antiferromagnet, while CeCoIns and Celrln5 
with smaller e than CeRhlns exhibit superconductivity. 
The difference in T c between CeCoIns and Celrlns may 
be attributed to the extent of two-dimensionality. Con- 
cerning the superconductivity and antiferromagnetism of 
CeTIns compounds, the combination of two dimensional- 
ity and the crystalline field splitting is necessary to reach 
a consistent picture. 
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FIG. 3. Phase diagram (a) in the U-U' plane for £=0 
and (b) in the U-e plane for U=U'. Meanings of ab- 
breviations used here are as follows: SC(F) denotes su- 
perconducting phase with P-symmetry, ICSO(q) indicates 
incommensurate spin ordered phase with wave vector q, 
ICOO(q) is incommensurate orbital ordered phase with 
wave vector q, and AFSO(Q) denotes antiferromagnetic 
spin ordered phase. Several wave vectors are defined 
as qi=(0.567r, 0.567r), q 2 =(0.537r, 0), q 3 =(7r,0), q x = (q x ,0), 
qy=(n,q y ), and Q=(7r,7r). Solid curves are the phase bound- 
aries determined by actual calculations, while dotted lines are 
schematic phase boundary guides for the eye. 

The orbitally degenerate model, Eq. (I), shows new, 
rich superconducting properties both for singlet and 
triplet pairings. In Fig. 3(a) and (b), we show the phase 
diagrams in the U-U' plane for e=0 and U-e plane for 
U=U', respectively. In these figures, four character- 
istic superconducting phases are observed: (1) a SF- 
mediated spin-singlet superconducting phase with Big- 
symmetry in the vicinity of a spin-ordered phase; (2) a 
g-wave superconducting phase with Aag-symmetry; (3) 
a spin-triplet superconducting phase with -E u -symmetry 
due to the cooperation between SF and OF around U=U' 
and e=0 and (4) an OF-mediated spin-singlet supercon- 
ducting phase with ^4i g -symmetry for U <CJ7' around an 
orbital-ordered phase. With respect to (1), it is thought 
that i?ig-superconductivity is induced by SF around 
q=Q=(7r,7r). Concerning the triplet superconductivity 
(3), note that for this region all superconducting insta- 
bilities are quite close to each other because of the multi- 
peak structures in % s (q) and x°(q) as shown in Fig. 1. 
We also note that the factors in front of the SF term in 
W s (q) are —3/2 for singlet and +1/2 for triplet pairing, 
while the factor for the OF term in IF°(q) is +1/2 for 
both pairings. Namely, OF is cooperative with SF for 



spin-triplet pairing, while SF and OF compete with each 
other for spin-singlet pairing 0. Thus, the spin-triplet 
pairing phase appears in the region for £«0 and U~U'. 
The mechanism of the OF-mcdiated singlet supercondu- 
tivity (4) is interesting. In this state, two quasi-particles 
with different pseudo-spin form an on-site intra-orbital 
singlet pair since U'^>U, leading to ^4i g superconductiv- 
ity. In short, singlet superconductivity with -Big- or A\ g - 
symmetry is stabilized when either SF or OF is dominant, 
while triplet superconductivity is favored when there is 
cooperation between SF and OF. 

In summary, we have studied superconductivity in the 
Ce-115 systems based on the orbitally degenerate Hub- 
bard model, and found that with increasing the tetrago- 
nal CEF splitting (4,2 _ y 2 -wave superconductivity due to 
AFSF emerges out of the suppression of OF in the vicin- 
ity of the AF phase. The concept of AFSF-induced su- 
perconductivity controlled by OF qualitatively explains 
the difference among Ce-115 materials. 
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